Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation

Background: Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth. Results: Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism. Conclusions: Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation

Cyanobacterial Cell Lineage Analysis of the Spatiotemporal hetR Expression Profile during Heterocyst Pattern Formation in Anabaena sp. PCC 7120

Diazotrophic heterocyst formation in the filamentous cyanobacterium, Anabaena sp. PCC 7120, is one of the simplest pattern formations known to occur in cell differentiation. Most previous studies on heterocyst patterning were based on statistical analysis using cells collected or observed at different times from a liquid culture, which would mask stochastic fluctuations affecting the process of pattern formation dynamics in a single bacterial filament. In order to analyze the spatiotemporal dynamics of heterocyst formation at the single filament level, here we developed a culture system to monitor simultaneously bacterial development, gene expression, and phycobilisome fluorescence. We also developed micro-liquid chamber arrays to analyze multiple Anabaena filaments at the same time. Cell lineage analyses demonstrated that the initial distributions of hetR::gfp and phycobilisome fluorescence signals at nitrogen step-down were not correlated with the resulting distribution of developed heterocysts. Time-lapse observations also revealed a dynamic hetR expression profile at the single-filament level, including transient upregulation accompanying cell division, which did not always lead to heterocyst development. In addition, some cells differentiated into heterocysts without cell division after nitrogen step-down, suggesting that cell division in the mother cells is not an essential requirement for heterocyst differentiation

Genomic Structure of an Economically Important Cyanobacterium, Arthrospira (Spirulina) platensis NIES-39

A filamentous non-N2-fixing cyanobacterium, Arthrospira (Spirulina) platensis, is an important organism for industrial applications and as a food supply. Almost the complete genome of A. platensis NIES-39 was determined in this study. The genome structure of A. platensis is estimated to be a single, circular chromosome of 6.8 Mb, based on optical mapping. Annotation of this 6.7 Mb sequence yielded 6630 protein-coding genes as well as two sets of rRNA genes and 40 tRNA genes. Of the protein-coding genes, 78% are similar to those of other organisms; the remaining 22% are currently unknown. A total 612 kb of the genome comprise group II introns, insertion sequences and some repetitive elements. Group I introns are located in a protein-coding region. Abundant restriction-modification systems were determined. Unique features in the gene composition were noted, particularly in a large number of genes for adenylate cyclase and haemolysin-like Ca2+-binding proteins and in chemotaxis proteins. Filament-specific genes were highlighted by comparative genomic analysis

Selection of Suitable Reference Genes for RT-qPCR Analyses in Cyanobacteria

Cyanobacteria are a group of photosynthetic prokaryotes that have a diverse morphology, minimal nutritional requirements and metabolic plasticity that has made them attractive organisms to use in biotechnological applications. The use of these organisms as cell factories requires the knowledge of their physiology and metabolism at a systems level. For the quantification of gene transcripts real-time quantitative polymerase chain reaction (RT-qPCR) is the standard technique. However, to obtain reliable RT-qPCR results the use and validation of reference genes is mandatory. Towards this goal we have selected and analyzed twelve candidate reference genes from three morphologically distinct cyanobacteria grown under routinely used laboratory conditions. The six genes exhibiting less variation in each organism were evaluated in terms of their expression stability using geNorm, NormFinder and BestKeeper. In addition, the minimum number of reference genes required for normalization was determined. Based on the three algorithms, we provide a list of genes for cyanobacterial RT-qPCR data normalization. To our knowledge, this is the first work on the validation of reference genes for cyanobacteria constituting a valuable starting point for future works

The Smallest Known Genomes of Multicellular and Toxic Cyanobacteria: Comparison, Minimal Gene Sets for Linked Traits and the Evolutionary Implications

Cyanobacterial morphology is diverse, ranging from unicellular spheres or rods to multicellular structures such as colonies and filaments. Multicellular species represent an evolutionary strategy to differentiate and compartmentalize certain metabolic functions for reproduction and nitrogen (N2) fixation into specialized cell types (e.g. akinetes, heterocysts and diazocytes). Only a few filamentous, differentiated cyanobacterial species, with genome sizes over 5 Mb, have been sequenced. We sequenced the genomes of two strains of closely related filamentous cyanobacterial species to yield further insights into the molecular basis of the traits of N2 fixation, filament formation and cell differentiation. Cylindrospermopsis raciborskii CS-505 is a cylindrospermopsin-producing strain from Australia, whereas Raphidiopsis brookii D9 from Brazil synthesizes neurotoxins associated with paralytic shellfish poisoning (PSP). Despite their different morphology, toxin composition and disjunct geographical distribution, these strains form a monophyletic group. With genome sizes of approximately 3.9 (CS-505) and 3.2 (D9) Mb, these are the smallest genomes described for free-living filamentous cyanobacteria. We observed remarkable gene order conservation (synteny) between these genomes despite the difference in repetitive element content, which accounts for most of the genome size difference between them. We show here that the strains share a specific set of 2539 genes with >90% average nucleotide identity. The fact that the CS-505 and D9 genomes are small and streamlined compared to those of other filamentous cyanobacterial species and the lack of the ability for heterocyst formation in strain D9 allowed us to define a core set of genes responsible for each trait in filamentous species. We presume that in strain D9 the ability to form proper heterocysts was secondarily lost together with N2 fixation capacity. Further comparisons to all available cyanobacterial genomes covering almost the entire evolutionary branch revealed a common minimal gene set for each of these cyanobacterial traits

The Transcriptional Landscape of the Photosynthetic Model Cyanobacterium Synechocystis sp. PCC6803.

Cyanobacteria exhibit a great capacity to adapt to different environmental conditions through changes in gene expression. Although this plasticity has been extensively studied in the model cyanobacterium Synechocystis sp. PCC 6803, a detailed analysis of the coordinated transcriptional adaption across varying conditions is lacking. Here, we report a meta-analysis of 756 individual microarray measurements conducted in 37 independent studies-the most comprehensive study of the Synechocystis transcriptome to date. Using stringent statistical evaluation, we characterized the coordinated adaptation of Synechocystis' gene expression on systems level. Evaluation of the data revealed that the photosynthetic apparatus is subjected to greater changes in expression than other cellular components. Nevertheless, network analyses indicated a significant degree of transcriptional coordination of photosynthesis and various metabolic processes, and revealed the tight co-regulation of components of photosystems I, II and phycobilisomes. Detailed inspection of the integrated data led to the discovery a variety of regulatory patterns and novel putative photosynthetic genes. Intriguingly, global clustering analyses suggested contrasting transcriptional response of metabolic and regulatory genes stress to conditions. The integrated Synechocystis transcriptome can be accessed and interactively analyzed via the CyanoEXpress website (http://cyanoexpress.sysbiolab.eu)